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We consider the effects of radiation damping on the electron dynamics in a Gaussian beam 
model of a laser field. For high intensities, i.e. with dimcnsionless intensity uq 3> 1, it is found 
that the dynamics divide into three regimes. For low energy electrons (low initial 7-factor, 
7o) the radiation damping effects are negligible. At higher energies, but still at 270 < ag, 
the damping alters the final displacement and the net energy change of the electron. For 
270 > flo one is in a regime of radiation reaction induced electron capture. This capture is 
found to be stable with respect to the spatial properties of the electron beam and results in 
a significant energy loss of the electrons. In this regime the plane wave model of the laser 
field provides a good description of the dynamics, whereas for lower energies the Gaussian 
beam and plane wave models differ significantly. Finally the dynamics are considered for 
the case of an XFEL field. It is found that the significantly lower intensities of such fields 
inhibits the damping effects. 
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I. INTRODUCTION 

The past few years have seen many advances in laser technology resulting in lasers of unprece- 
dented powers and intensities, the current record being approximately lO^^W/cm^ [Ij. This trend 
is expected to continue throughout the next few years, such as with the European Extreme Light 
Infrastructure (ELI), a facility that may deliver intensities as high as lO^^W/cm^ [2]. Such extreme 
intensities will allow the probing of fundamental physics in previously inaccessible regimes. 

In fields of such high intensities, an electron will be accelerated so strongly that its own radiation 
emission may significantly affect its motion. This opens up the possibility of testing experimentally 
the classical theory of radiation reaction in such a context. Moreover, even when the classical theory 
is no longer directly relevant, such as when considering laser induced strong field QED processes, 
an understanding of radiation damping effects is still important. This is because many such QED 
processes either have close classical analogues (e.g. nonlinear Compton scattering pT, ^3]) or the 
current models describing them incorporate classical theory. An example of the latter is the topical 
area of runaway QED cascading in strong laser pulses. The numerical modelling of this phenomena 
often involves propagating the particles classically through the laser, whilst altering the particle 
number at each time step via quantum transition rates (see e.g. OEj)- 

An important factor that needs to be taken into consideration is the description we choose 
for the laser pulse. Many recent discussions of radiation damping effects have restricted their 
consideration to plane wave models of the laser field [JHS]. However, the new generation of high 
intensity facilities will, in part, achieve their high outputs by a strong focussing of the laser pulse. 
This suggests that the spatial properties of the beam will become more important and so we should 
adopt a more realistic description of the field. In this paper will will use a Gaussian paraxial 
beam approximation, including terms up to fifth order in the expansion parameter. This beam 
model describes the spatial as well as the temporal properties of the beam, giving us an accurate 
description of the field profile. We will also compare our results with those of the plane wave model, 
which will allow us to evaluate the accuracy of such a description. 

In this paper we have two aims. Firstly, to assess the significance of radiation damping on the 
electron dynamics in a realistic beam profile. Secondly, to look for regimes where the radiation 
reaction effects are particularly prominent, allowing for an easy experiment testing of the classical 
theory. 
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II. THEORY 
A. Equation of motion 

A charged particle undergoing acceleration (e.g. due to the presence of a laser field) will emit 
electromagnetic radiation. The emission of this radiation will result in a momentum loss for the 
particle, thus subsequently affecting its motion. The classical action describing such a system is 

mm 

S = -mjdT-ej (fxfA^ ~ \ j xF^uF^\ (1) 

where m is the mass (from hereon taken to be the electron mass), e the charge, the four-current 
and we have adopted units where the speed of light is unity, c = 1. The gauge potential and 
electromagnetic field tensor are denoted by A^^ and F^'^ respectively. Varying ([T]) with respect to 
and gives us the governing equations 

d^F>^^ = f, (2) 
mu^' = eF^"'uy = F^^. (3) 

In the second equation is the particle four-velocity and the dot denotes differentiation with 
respect to proper time r. The electromagnetic tensor F^^^ may be decomposed into the sum of the 
external field (in this case the laser field) Fji^^ and the radiation field emitted by the electron F'^^ 

F^^ = FZ + FZ- (4) 

The resulting equation of motion, after calculating F^^ is the Lorentz Abraham Dirac (LAD) 
equation [T2HIi] 

2 

mil'' = eFZuu -~ iuf'u'' - u'^il'') u^. (5) 

ext 3^^v V / 

This equation is infamous due to the presence of the u terms on the right hand side which give 
rise to (unphysical) runaway solutions. A common way to remove this problem is to approximate 
u using the first term in ^ (i.e. the Lorentz force expression). This results in the Landau Lifshitz 
(LL) equation [llj 

f 2 2 

= ^F'^'^u, + ro <^ ^F'^'^u, + ^F'^^F^'^u, - ^u^F'^^Fjup 

which is a perturbative expansion of LAD to first order in the coupling rg = 2a/3, where a = e^/47r 
is the fine structure constant, and has the benefit of no longer exhibiting the runaway solutions. 
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For a laser beam described by the wave vector we define a dimensionless measure of laser 
intensity (see [T5j ) 

"■0 ~ 2 n \2 ' vV 
[k ■ pqY 

where pg = ttiUq is the initial electron four-momentum and the brackets (...) denote the maximum 
value of the enclosed quantity over a laser cycle. For a plane wave laser in the lab frame, this takes 
the more familiar form 

ao = , (8) 

ujm 

where oj is the laser frequency and Aq the magnitude of the electrical field strength. 

To neglect the influence of radiation damping we set tq = 0. Then the equation of motion ([g]) 
takes the form 

mu^' = eF^^'^Uy, (9) 

which is, of course, the Lorentz force equation. Using the solution to ([9]), the standard measure of 
radiated energy loss is given by Larmor's formulsj^ 

P = rov?, (10) 

where P is the radiated power. The energy loss of the electron over a laser cycle is given by 

R= — =ro-al^{l + p). (11) 
bjm m 

When this parameter reaches unity we are in the "radiation dominated regime" |17j where the 
radiation damping effects are of the same magnitude as the Lorentz force. In order for the LL 
equation to be valid as a perturbative expansion of LAD we require that the radiation damping 
term is much smaller than the Lorentz force term. Since it is the final term of ([g]) that gives us the 
biggest contribution to the damping, we require that 

roe^UaF'^^'F^ ^up /m^ < eF^^'uylm, (12) 

which, for a head-on collision in the lab frame between the electron and a plane wave laser field, 
gives us the constraint 

rou;^ao7^/m ^ a;ao7. (13) 



^ Using just the Lorentz force solution for u is consistent with the LL equation being an expansion truncated to first 
order in ro |16| . 
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For an optical laser this becomes 

ao7^ < 10^ (14) 

Indeed, in [S] the authors find that the solutions to the LAD and LL expressions begin to differ 
even at an order of magnitude below this estimate. 

As well as ensuring that we stay in a regime where the derivation of the LL equation is valid, 
we must also take a moment to ensure that we keep within the domain of classical physics. A 
common measure of the importance of quantum effects is to consider the work done by the laser 
field over the distance of a Compton wavelength 

X = 3 • (15) 

This was first introduced in jl8| [T9] and is expressed here in units where h = 1. Quantum effects 
become important when this measure approaches unity. In the lab frame we find that to be in the 
classical regime we require 

ao^oj ^ m, (16) 



which, for an optical laser, can be a more stringent requirement than (14) 



B. Description of the field 

We describe our laser field using a pulsed Gaussian beam model, in a similar way to Salamin et 
al in [20j. Proceeding in the same way as McDonald |2T] we begin by adopting the Lorentz gauge 

and dictate that any vector potential describing our laser field must satisfy the vacuum wave 
equation 

V^.^l^. (IS, 

We take the laser to propagate in the +z direction and assume a generic potential, linearly polarised 
in X 

A = xAogir])i;ix,y,z)e-''\ (19) 
where Aq is the wave amplitude, rj = ujt — kz, and is a generic pulse shape function. Inserting 



( 19 ) into ( 18 ) gives us 



V^V; - 2ifc|^ (^1 - / ) = 0, (20) 
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where g' = dg/drj. In general it is hard to satisfy (20) since ■0 is a function of (x,y,z) and g is a 



function of the phase ry. To proceed we begin by rescaling our coordinates 

^ X y ^ z 

C = — , = — , ( = —, 

Wq Wo Zr 



(21) 



making them dimensionless. Here wq is the beam waist diameter and Zr = kwQ/2 is the Rayleigh 
length. Following [21] we specify that the pulse shape function satisfies 



9 <-g- 



Equation ( 20 ) can then be approximated by 



Vi^-4^.o^g^O, 



(22) 



(23) 



= — + — , ^ = 7/;(e,l/,C), 



(24) 



where 

q2 pa 

and we have introduced the aspect ratio = wo/zr which, when small, closely approximates the 
beam diffraction angle. Since is typically small, we can expand ip in the series 



i) = i)Q + 9111)2 + ^^04 + . . 



Equating coefficients of we have, from (23) 



etc. 



(25) 

(26) 
(27) 
(28) 



Equation (26) is the paraxial wave equation from traditional Gaussian beam theory. Its solution 



is the well known first order Gaussian beam solution 
where 

/ 



1 



giarctanC^ pi = f ^ ^2 _ 



The solution to (27) was originally found by Davis [22] 

7 , fp' 



i'2 



2+ 4 



V'o, 



(29) 



(30) 



(31) 



and Barton and Alexander [23] proceeded to find the solution to ( 28 ) 



(32) 



Analogously to the case of the vector potential (19), we assume that the scalar potential can be 
written in the form 



(p{t, X, y, z) = g{vi)^{x, y, z)e*''. 



(33) 



Then the Lorentz gauge condition (17) gives us 



which means that 



„ iu}(j) 1 — i— 1 iuj<i 

dt V g 



-V- A. 



(34) 



(35) 



Thus our electric and magnetic field components may be found from (19) via 



E = -ikA- -V(V- A) 



k 



B = V X A, 



(36) 
(37) 



(for details of the calculation see |20^ I23j). Taking the real part of the resulting expressions gives 
us (to fifth order in Oq) 

+i [4*52 - Sp^Ss - 2p2(p2 _ 16^2)54 - 4/(/>2 + 2^2)^^ + ^8^^] ^ ^ (3g) 
Ey = Piu {elS2 + I [4/52^4 - /Ss] ) , (39) 



E^ = Pi( 9oCi + I [-2C2 + 4p2c73 - /C4] 



+ ^ [-I2C3 - 12p2C4 + 3VC5 - 12p6C6 + p^Cy] 



5, = 0, 



(40) 
(41) 



By = p(^So + ^ [2p^S2 - p'Ss] + ^ [-4S2 + Sp^Ss + Wp^S4 - Sp^Ss + p^Se] ) , (42) 
= Pu^e^Ci + ^ [2C2 + 2p2C3 - 

(43) 



+ ^ [12C3 + 12^2^4 + 6/C5 - Sp^Ce + p^Ct\ ] , 
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where the prefactor is given by 



P = Ao— 5(r?)exp(-4), = + y'' . (44) 



Here w = vj{z) \s a. measure of the beam diameter according to 



, 2 

z 



w{z)=wo\l+ - , (45) 



and the functions Sj and Cj are defined 



5j = (^ysine, (46) 



w I 



Cj = (^ycose, (47) 



w J 



where 



e = r/ - — + (i + 1) arctan C, (48) 

where H = z + z"^/ z \s the radius of curvature of the field. 

It has been shown in [20j that, to accurately describe the properties of an optical beam for 
realistic parameter values {wq ~ 5/um), it is necessary to keep terms up to at least 0{9'q). With 
this in mind we will henceforth work with the full expressions (to order 6^) given above. (However, 
in the case of an X-ray free electron laser (XFEL) the beam is much less tightly focused, typically 
having a waist that is many wavelengths in diameter. In such cases the expansion parameter, 9q, 
is small enough that one only needs to use the first term without corrections.) All that remains is 
to define our pulse shape function g{r]). A common choice when studying plane wave dynamics is 
to use a Gaussian pulse profile g = exp[— (r^/r/o)^], where r/o is the number of cycles in the pulse. 



However, this function does not always satisfy the constraint (22). Instead, in this study we adopt 
the profile recommended by McDonald 



5(77) = sech^^^. (49) 
In Figure [1] we show the laser pulse intensity for typical parameter values and at various times. 



C. The set-up 

Our aim is both to assess the significance of radiation damping in different regimes and to focus 
on those regimes where the effects are most prominent. The latter cases will allow for experimental 
tests of the classical theory. We see from ^ that the coupling tq between the radiation damping 
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FIG. 1. Plot showing the energy density {{E^ + B'^)/2) at different times. A: ujt = -2000, B: cut = -1000, 
C: Ljt = +0.5, D: ut = +1000, E: ut = +2000. Parameters are A = l^m, wq S^m, % = 20. The sohd 
black lines show the boundary of the beam waist wq. 



force and the Lorentz force is very small. In order to compensate for this we need either for the 
damping force to be large, or for the electron to experience the damping force for a long time 
period. For the damping force to be large the electron must undergo a large acceleration. In 
practice this means that the electron must see a laser field of high intensity - either by having a 
large oq or by having a head on (or near head on) collision between a moderately high intensity laser 
and a moderately high energy electron (moderately high meaning that we adhere to the physical 
constraints ( 14|16| )). The other option is for the electron to spend a longer time period in the 
laser field. This is best achieved via a capture and acceleration scenario (CAS) |24H26j where the 
electron enters the field and gets captured near the focus and subsequently accelerated. For such a 
situation to occur we should have a high energy electron interacting with the laser pulse in a same 
direction (or near same direction) collision. Although in the electron's frame the laser is not as 



10 



strong as it would be for a head on collision, the fact that the electron is fast moving means that 
it will travel along with the laser pulse for a significant period of time and so, even though they 
are smaller, the damping effects will accumulate. The initial conditions must be chosen such that 
the electron is captured by the beam, rather than being ejected due to the ponderomotive force. 
This is generally achieved by having a suitable balance between oq and the initial electron energy. 
It is also advantageous to have our pulse length r]Q as large as possible (but, of course, the pulse 
length is constrained by the total amount of energy the laser facilities are able to deliver). 

In Figure[2]we give a schematic of our set up. The simulation will take place in three dimensions, 
but for simplicity we will consider things here in the y = plane. We start our simulation when 
the electron is at the point {zq,xq) with energy 70 moving towards the point (zi,0), making an 
angle 6 with the beam axis. When zi = the electron is aimed at the centre of the pulse, otherwise 
{zi 7^ 0) the electron is sent towards the right or left of the focus. The beam is propagating in the 
+z direction; therefore collisions where ~ will be referred to as 'near same direction' and cases 
where 9 ~ 180 'near head on'. 

We will solve the equation of motion ^ numerically using a second order leapfrog method. To 
make the process less computationally expensive, we exclude the first of the radiative correction 
terms (the term proportional to F^'^) from our calculation. This term goes like roao7 which, for an 



optical laser, even at the maximally allowed value of ao7 (16) is still 2 orders of magnitude smaller 
than the Lorentz force term. Test runs of our code indicate that this is a justified approximation. 



III. RESULTS 
A. Plane wave dynamics 

Before considering the full Gaussian beam model, we will first devote our attention to a simpler 
case ~ that of the plane wave. We obtain plane wave fields in the limit where the beam waist 
becomes large, wq ^ 00. Then the laser field takes the configuration 

E^ = Aog{v)smr], Ey = 0, E, = 0, (50) 
= 0,By = Aog{r]) sin i], = 0. (51) 

Plane wave models are commonly used when modelling laser-electron interactions since they are 
simple enough to allow the properties of many relevant physical phenomena to be calculated ana- 
lytically, whilst still retaining important characteristics of the laser field (such as time-dependence). 
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FIG. 2. Schematic of our simulation set up shown in the y — plane. The laser is propagating in the 
k = +z direction and the electron starts at (zq^xq) heading towards (zi^O) at an angle 9 to the z axis. 

From our perspective, an additional reason for considering plane waves is that the analytical solu- 
tion to the LL equation is known |27j . thus allowing us to benchmark our codes. However, when 
working with plane wave models one must proceed with caution. Plane waves are infinite in their 
spatial extent and therefore unsuited to situations where the laser-electron interaction occurs over 
a long spatial distance. For example, a high energy electron propagating in the same direction as 
the laser will spend a long time in the laser pulse since it will be travelling at nearly the speed of 
light. In such circumstances a long interaction time equates to a long spatial distance of interaction 
in the propagation direction. It can be seen from 
z-direction, a realistic laser field will become significantly damped (decaying like 1/w ~ 1 — z'^/z^), 
whereas the plane wave model will maintain a constant peak amplitude throughout. Because of 
this we will only consider the case of a near head on collision in our plane wave analysis. 

The left hand plots of Figure show the case of a near head on collision [6 = 170°) between an 



38||43) that, even over a modest distance in the 
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FIG. 3. Near head on collision {9 — 170 degrees) for plane wave laser field. Left: gq = 150, 70 = 100, right: 
ao — 250, 7o = 150; 770 — 20 for both cases. Top plots show trajectories and bottom plots show corresponding 
electron energies (7-factors). Black lines: undamped solutions (Lorentz force), red lines: damped solutions 



electron of initial energy 70 = 100 and a plane wave of intensity oq = 150. The pulse parameter is 
chosen such that the laser has 20 cycles (770 = 20) which, for the 1/im wavelength we are considering, 
corresponds to a pulse duration of approximately lOfs. We can see that for these parameter values 
radiation damping only has a minimal effect on the electron's trajectory. Nevertheless, as can be 
seen from the lower left plot, the electron loses approximately 50% of its energy due to radiation 
damping. This is consistent with the analysis presented by Koga et al |28| (see also [9]). The right 
hand plots of Figure |3] show the same quantities but for different parameter values - this time 
ao = 250, 70 = 150. These parameters are towards the higher end of the range allowed by our 



conditions (14, 16). What is most noticeable in these plots is that the radiation damping is now 
strong enough to reflect the electron so that for a while it co-propagates with the laser. The onset 
of this reflection regime was first identified by Di Piazza et al as a means to experimentally 
test the LL equation via the resulting Thomson scattering spectra. For a head on collision (i.e. 
9 = 180°) between an electron and a short plane wave laser pulse the authors derive an expression 
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FIG. 4. Left plots: longitudinal displacement D of the electron as a function of 70 for a plane wave, near 
same direction collision {0 — 170 degrees). Right plots: net energy change — 70 as a function of 70. Top 
plots have parameters ao = 180, 770 — 20 and bottom plots oq ~ 250, 770 = 20. Black lines: undamped 
solutions, red lines: damped solutions, blue dotted lines indicate 1% of i?, blue dashed lines: 270 = oq. 



for the onset of reflection. This is 



R> 



2al 



> 0, ao > 1, 



(52) 



which, below the radiation dominated regime R, reduces to 270 > aojf] This is broadly consistent 

with our results, bearing in mind that we are not considering a directly head on collision. 

To investigate this phenomena more fully we will fix our intensities at oq = 180 and 250 and (for 

both the undamped and the damped solutions) determine numerically the size of the longitudinal 

shift of the electron as a function of 70. We define our longitudinal displacement parameter D 

^ Further light may be shed on this by considering the kinematics in terms of a laser induced mass shift of the 
electron. If one averages the electron momentum over a laser cycle, one finds that the electron can be considered 
to have undergone a mass shift rn^ — )■ m? {1 + al). Then one finds that when qq ~ 270 the laboratory frame can 
be interpreted as an intensity-dependent centre-of-mass frame. See [3^ for a fuller discussion. 
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to be the longitudinal shift that the electron undergoes before it reaches a distance |x| = 20/um 
from the laser axis, compared to where it would have been if no field was present. (It should 
be noted that our displacement parameter D measures a quantity that is not directly analogous 
to the 'reflection' referred to in |17] . This is because our D parameter measures the fleld induced 
change in longitudinal position of the electron at a distance from the beam axis. Thus this measure 
registers angular deflections as well as direct reflections, which is precisely what a detector would 
measure if it were placed alongside the beam.) The displacement parameter is shown in the left 
hand plots of Figure [4] and the net energy gain (7/ — 70) on the right. 

At low energies the electrons are strongly displaced/reflected, but experience only a minimal 
change in energy. A useful criteria for demarking this region is found to be 70 such that R = 0.01, 
i.e. 1% of the radiation dominated regime. As we increase 70 above this value we flnd that the 
radiation damping begins to play a role in the dynamics, with the damped electron being subject 
to a signiflcantly higher displacement than its undamped counterpart. In this regime the damped 
electron loses a large proportion of its initial energy, while the undamped one remains unaffected. 
Beyond this, in the reflection regime 270 > ao, it is only the damped electron that experiences a 
displacement. This, together with our results in Figure [3j indicate that we have entered a regime 
of radiation reaction induced electron capture. 

B. Gaussian beam dynamics 

Having throughly considered the plane wave dynamics, we now turn our attention to the more 
realistic pulsed Gaussian beam model of the laser. We begin by looking at the case of a near same 
direction collision. As stated earlier, in this case an electron in its rest frame will see a laser of 
lower intensity than if it were counter propagating with the beam. The advantage of this set up is 
that it allows for a CAS scenario where the electron will get captured and carried along with the 
pulse. Thus radiation damping effects, although somewhat smaller, may accumulate over a longer 
time period. Choosing our parameters carefully so that CAS takes place, one finds that, even 
though the energy gains of the electron are large (see Figure [S]) , the effect of radiation damping is 
negligible. This is true even for the very highest parameter values allowed in the classical regime 
and is consistent with the results of Mao et al [29] 

Therefore we now turn our focus to the case of a near head on collision. In Figure [6] we show the 
electron trajectories and energies for the same parameter values as the plane wave case in Figure 
Isl For the more modest parameter values in the left hand plots (oq = 150, 70 = 100) we find 
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FIG. 5. Near same direction collision [9 — 10 degrees) between an electron and a pulsed Gaussian beam. 
Parameters are ao — 600, 7 ~ 400, 770 = 20. Black lines (solid): undamped solutions, red lines: damped 
solutions. Black lines (dashed): beam waist. 

that the behaviour of both the damped and undamped electrons is very similar to the plane wave 
case, with little reflection and radiation damping causing an energy loss of roughly 50%. At the 
higher parameter values in the right hand plots (ag = 250, 70 = 150 - which are again towards the 
higher end of what is allowed in the classical regime) the dynamics have altered somewhat from the 
plane wave case. Although the damped electron is reflected by the beam, it leaves it in a different 
direction. Nevertheless the electron 7-factor has a very similar behaviour to the plane wave case. 

As with our plane wave analysis, it is instructive to fix the laser intensity and consider the 
longitudinal displacement and final electron energy as a function of 70. The results of this analysis 
are shown in Figure [7} Once again there are three distinct regimes, although the criteria we used 
to demarcate the first of these regions in the plane wave case is now only suitable as a rough guide. 
Generally, we find that the radiation damping effects kick in later than for the plane wave case. This 
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FIG. 6. Near head on collision {9 — 170 degrees). Left plots; uq — 150, 70 — 100; right plots: gq — 250, 
7o — 150; 770 — 20 for both cases. Top plots show trajectories and bottom plots show corresponding electron 
energies. Black lines (solid): undamped solutions, red lines: damped solutions, black lines (dashed): beam 
waist w. 



can be seen from the plots by observing tliat the damped and undamped lines are still together 
beyond the l%i? limit (which is where they separated in our plane wave examples). That the 
radiation damping effects occur later for the Gaussian pulse field is to be expected, since the field 
intensity decays more quickly away from the focal point than in the plane wave model. Observe 
that in the case of low 70 both the damped and undamped electrons experience a net energy gain, 
whereas in the plane wave case it was only the damped electron that underwent an energy change. 
The marker 270 = oq still serves as a good indicator as to when the undamped electron can pass 
through the beam without being displaced. Therefore, once again, in the regime 270 > qq we have 
a situation where the damped electron is displaced but the undamped one is not. Examining the 
trajectories (e.g. Figure |6| we see that we are in a regime of radiation reaction induced electron 
capture. Comparing Figures [7] and |4] indicates that in this regime the displacements and energy 
gains are well described by the plane wave model. To investigate further, in Figure [8] we plot the 
relative error of the plane wave approximation to the final energy (7Gauss — 7Pw)/7Gauss- Beyond 
270 = ao the plane wave model gives us an extremely good measure of the net energy change, even 
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FIG. 7. Left plots: longitudinal displacement D of the electron as a function of 70 for a Gaussian pulse, 
near same direction collision {6 ~ 170 degrees). Right plots: Net energy change 7/ — 70 as a function of 70. 
Top plots have parameters ao — 180, 770 = 20 and bottom plots ao — 250, 770 = 20 . Black lines: undamped 
solutions, red lines: damped solutions, blue dotted lines indicate 1% of i?, blue dashed lines: 270 = oq. 



though the actual dynamics are a httle different. That the two models converge can be understood 
by recalhng that the high intensity limit is the same as the low frequency limit. In the low frequency 
limit both the plane wave and Gaussian beam models behave as constant crossed fields. 

In a realistic situation it will not be possible to ensure that the electrons will all be aiming 
precisely for the centre of the beam focus. For our results to be interesting experimentally, we 
need to consider their sensitivity to the spatial size of the electron beam. Therefore we conduct 
some simulations where small bunches of electrons interact with the beam. A typical source of 
electrons in experiments such as these is from a linac. In our simulations we will model the electron 
beam as being 10/um in diameter, which is typical of facilities such as the ELBE accelerator at 
the Forschungszentrum Dresden- Rossendorf in Germany ^30j . We will assume the beam to be of a 
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FIG. 8. Relative error of the final electron energy 7/ calculated using the plane wave approximation and 
the full Gaussian beam expression, (7Gauss — 7Pw)/7Gauss- Parameters are oq = 180, rjo = 20. Black lines: 
undamped solution, red lines: damped solution, blue dotted line: 270 = uq. 




FIG. 9. Near head on collision {6 = 170 degrees). Parameters are Oq = 250, 7 = 100, rjQ = 20. Black lines 
(solid): undamped solutions, red lines: damped solutions, black lines (dashed): beam waist w. 

high quality with negligible energy spread (e.g. A7/70 = 10^^ is feasible with the aforementioned 
facility [3]) and that the concentration of electrons is dilute enough for Coulomb repulsion between 
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the electrons to be neglected. The bunches will be centered at the beam focus, extending to 5/um 
either side of it (i.e. zi = — 5/im. . . 5/im), with all the electrons travelling parallel to one another. 
The resulting particle trajectories with and without damping are shown in Figures [9]and[T0| In 
Figure [9] we choose our parameter values (70 = 100, oq = 250) from the region [R > 1%, 270 < ao]) 
i.e. the regime where both the damped and undamped electrons experience displacement (see 
Figure [7]). One sees that the damped and undamped dynamics are very different to one another. 
The scattering direction of the undamped electrons is much more sensitive to zi than for the 
damped electrons. The behaviour of the damped electrons is more consistent, with all of them 
being reflected, although the exit direction flips as we move slightly in front of the laser focus. 



For the second set of parameter values (Figure 10, 70 = 150, oq = 250) we are in the radiation 



reaction induced capture regime (270 > ao) where the electrons have sufficient energy so that the 
undamped solutions exhibit only minimal displacement, but radiation damping is strong enough 
to capture the damped electrons. In this regime the undamped electrons are well culminated, 
exhibiting much less deflection than in the previous example and, once again, all of the damped 
electrons are captured by the beam. We find that, although there is still a flip in the direction of 
the ejected electrons, the essential phenomena of radiation reaction induced capture is stable with 
respect to electron beam width effects. The net energy change of the damped electrons is a little 
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FIG. 11. Final gamma as a function of zi. Near head on collision {9 — 170 degrees). Parameters are 
flo = 250, 7 = 150, 770 = 20. Black lines/circles: undamped solutions, red lines/squares: damped solutions. 

more sensitive to zi, however, as can be seen by considering the final 7- factors of these electrons 



(Figure 11). 



C. XFELs 



Let us now briefly consider how the situation changes if we swap our optical laser fleld for an 
XFEL beam. The properties of XFEL beams are somewhat different to the optical fields we have 
been considering [31j, and this will reflect itself in the subsequent dynamics. As well as the obvious 
shorter wavelength, the main differences are the weaker focusing of the beam and the much higher 
number of cycles in the pulse. The weaker focusing means that the beam waist is many laser 
wavelengths in diameter. A typical lOOfs pulse with wavelength A = Inm will contain of order 10^ 
cycles. Thus the XFEL pulse behaves almost as if it were a continuous beam with stationary focus 
(i.e. 5 ~ 1). This means that the electron will spend a longer time in a weaker part of the pulse 
before reaching the focus, thus needing a comparatively higher initial energy in order to probe the 
most intense part of the beam. The upshot of this is that the spatial effects of the beam will play 
a bigger role, which may mask some of the radiation damping effects. 
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In Figure 12 we show the results of a numerical study of the dynamics of electrons in typical 
XFEL fields. The left hand plots show the electron displacement and net energy change for the 
case of a Inm wavelength field of intensity oq = 0.002. This intensity corresponds to lO^^W/cm^ 
and is likely to be the highest that can be achieved with the newest generation of facilities, such 
as the European XFEL [31J. It is clear from the plots that, even though we are far into the 
regime 270 > oq, the intensity is just not high enough for radiation damping to be significant. 
We could of course increase the electron energy but, since the electron energy would then be far 
greater than the laser intensity, the interaction time would be very short. Also, and perhaps more 
importantly, the derivation of the LL equation is not valid for higher 70 with these parameter 
values (see (|14[)) and so we would have to find an alternative approach to describing the dynamics. 
Instead we consider the dynamics with a larger oq . The right hand plots of Figure [12] show the 
same quantities but with uq = 2. For these values the effects of radiation damping are once again 
significant, resulting in an increased displacement and net energy loss to the electrons. However, 
upon examining the trajectories we find that the difference in displacement is caused by a deflection 
of the damped electron rather than a reflection. Therefore we don't have the same phenomena of 
radiation reaction induced electron capture as we did in the optical case. 



IV. CONCLUSIONS AND SUMMARY 



We have analysed the impact of radiation damping on the dynamics of electrons in high intensity 
pulsed laser fields. We began by considering an idealised pulsed plane wave model of the laser field 
in a 'near head on' (9 = 170 degrees) collision with a counter propagating electron. For relatively 
low parameter values (70 = 100, ao = 150) we confirmed that the radiation damping causes an 
energy loss to the electron (in this case ~ 50%) but has only a minimal impact on its trajectory. At 
high intensities it was found that the damped electron can become refiected by the beam and we 
saw evidence of radiation damping induced electron capture. Performing a more detailed analysis, 
we found that the dynamics fall neatly into three regimes. For low 70, such that R < 1%, the 
electrons are strongly deflected by the pulse but the effects of radiation damping are negligible. 
If we increase 70 such that we are in the range [R > l%,27o < ao], then damping effects come 
into play. Both damped and undamped electrons undergo displacement, but the damped electron 
experiences a net energy loss. In the very high energy regime, 270 > oq, the damped electron 
is displaced but the undamped one is not. We then moved to the more realistic paraxial pulsed 
Gaussian beam model of the laser field. This allowed us to take into account the spatial as well 
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FIG. 12. Study of the dynamics in an XFEL beam. Parameters are wq — 0.05/im, A — Inm, 6 = 170 
degrees and rj^ = 1.89 x 10''', implying a pulse duration of ^ lOOfs. Left plots: ag — 0.002, right plots: 
flo — 2. Top plots: longitudinal displacement D of the electron as a function of 70 (Note that the values are 
calculated for a point 2/im from the beam axis, rather than 20/im as they were in the optical examples.). 
Bottom plots: net energy change 7/ — 70 as a function of 70. Black lines: undamped solutions, red lines: 
damped solutions. 

as the temporal effects of the beam. Taking our beam to be focused to a waist of wq = 5/.im, we 
considered two different intensities - ao = 180 and 250 - and studied how the subsequent dynamics 
change as a function of the initial electron energy 70. We found once again that the dynamics 
fall into three regimes. For low 70 the radiation damping effects are again negligible. However, 
damping effects start to become important at higher energies than for the plane wave examples, 
with R = 1% no longer being a reliable indicator of the onset of the damping regime. Also unlike 
the plane wave case, for low 70 both the damped and undamped electrons in the Gaussian pulse 
experience a net energy gain. In the intermediate region [70 ^ 1,270 < ao] both the damped and 
undamped electrons were found to be displaced/deflected by the beam. The condition 270 = ao 
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once again provided a good criteria for when the undamped electron is no longer deflected. In the 
high energy regime 270 > ao it was found that only the damped electron gets displaced. Thus one 
can identify this regime as one of radiation damping induced electron capture. Another feature 
of this regime is that the electron dynamics are qualitatively the same for the Gaussian beam as 
for the plane wave, with the plane wave model providing an extremely good approximation to the 
net energy change. Also, giving consideration to the spatial effects of the electron beam, we found 
the dynamics to be reasonably stable with respect to changes in the electron position. Finally, 
radiation damping in XFEL type fields was briefly considered. It was found that the relatively low 
intensities expected in the current generation of facilities will limit any damping effects. 

Thus we have analysed the effects of radiation damping in a realistic model of a laser field. What 
remains is for our results to be tested experimentally. An experimental test of our radiation reaction 
induced electron capture scenario, in particular, will allow an evaluation of the applicability of the 
classical theory in the context of such high intensity laser fields. 
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